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Abstract. Recently, water absorbed in the porous silica material MCM-41-S15 has 
been used to demonstrate an apparent fragile to strong dynamical cross-over on cooling 
below '-^220K, and also to claim that the density of confined water reaches a minimum 
at a temperature around 200K. Both of these behaviours are purported to arise from the 
crossing of a Widom line above a conjectured liquid-liquid critical point in bulk water. 
Here it is shown that traditional estimates of the pore diameter in this porous silica 
material (of order 15A) are too small to allow the amount of water that is absorbed 
'^ I by these materials (around 0.5gH2O/g substrate) to occur only inside the pore. Either 

the additional water is absorbed on the surface of the silica particles and outside the 
Q ', pores, or else the pores are larger than the traditional estimates. In addition the low 

Q Bragg intensities from a sample of MCM-41-S15 porous silica under different dry 
^SJ ■ and wet conditions and with different hydrogen isotopes are simulated using a simple 

^ I model of the water and silica density profile across the pore. It is found the best 

agreement of these intensities with experimental data is shown by assuming the much 
tH- . larger pore diameter of 25A(radius 12. SA). Qualitative agreement is found between 

CO ' these simulated density profiles and those found in recent empirical potential structure 

f~^ ■ refinement (EPSR) simulations of the same data, even though the latter analysis did 

^3 I not specifically include the Bragg peaks in the atomistic structure refinement. It is 

shown that the change in the (100) peak intensity on cooling from 300K to 210K, which 

previously has been ascribed to a change in density of the confined water on cooling, 

can equally be ascribed to a change in density profile at constant average density. It is 

j^ ' further pointed out that, independent of whether the pore diameter really is as large 

as 25A or whether a significant amount of water is absorbed outside the pore, the 
earlier reports of a dynamic cross-over in supercooled confined water could in fact be 
a crystallisation transition in the larger pore or surface water. 
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1. Introduction 

In a series of papers which extend back at least 8 years, S-H Chen and colleagues 
have performed detailed studies of the structural and dynamical properties of confined 

water . p El [31 II El El El [91 [ini ni la [131 [11 [la 23 El [Tg [ig ^ 

include quasi-elastic neutron scattering (QENS), small angle scattering (SANS), NMR, 
and computer simulation, and there is a general finding that confined water appears to 
undergo an dynamical crossover from Arrhenius to Vogel-Fulcher (fragile to strong, FS) 
behaviour at a temperature near 220K. The hypothesis is made that this FS transition 
coincides with the crossing of a Widom line in the phase diagram of supercooled water 
which extends above a purported liquid-liquid critical point [3j . The same hypothesis is 
made concerning the apparent observation of a density minimum in confined water at 
a lower temperature near 200K |22], with the lowest density occurring at temperatures 
below that of the Widom line, and the maximum decrease in density occurring at the 
Widom line itself. However the dynamic crossover phenomenon is not quite a universal 
phenomenon of confined water since in some cases, where the surface is hydrophobic, 
the effect appears to be absent or shifted to a much lower temperature regime. p3l [24] 
suggesting it is impacted to some extent by the nature of the surface itself. 

Others have questioned the interpretation of the dynamic data in terms of a dynamic 
crossover. Both Cerveny et. al. [25] and Swenson [26] issued comments on Liu et. 
a/.[3j. Cerveny et. al. argue that the observed fragile-to-strong crossover is in fact due 
to the onset of confinement effects and quote several related cases of confined water 
where the same behaviour in the QENS data has been seen, and also quote the case of 
polymer blends where the same trend is seen. Swenson argues that if the QENS data are 
taken literally they would extrapolate to a glass transition temperature of 50K, which 
would be unacceptably low. Like Cerveny et. al. he also argues that what the QENS 
data are seeing is the effect of confinement killing the a relaxation process in water, 
rather than any FS transition. Subsequently Swenson and coworkers have published a 
dielectric relaxation of study of water highly confined in MCM-41 [27]. They observe no 
obvious transition from Arrhenius to Vogel-Fulcher- Tammann (VFT) behaviour in the 
dielectric relaxation time of the water as a function of temperature, particularly at the 
temperature where this transition occurs in the QENS and NMR data. In a more recent 
study of water near a protein surface Doster et al. also suggest an alternative explanation 
based on a glass transition scenario, instead of the purported fragile-strong crossover. [28] 
Even more recently Limmer and Chandler have questioned whether atomistic computer 
simulations of low temperature water are capable of observing a liquid-liquid transition, 
arguing instead that what is observed in the simulations is a non-equilibriated liquid- 
crystal transition [29] . 

To try to understand what might be going on with the structure of the water 
confined in a pore of MCM-41-S15 Mancinelli et al. undertook a combined neutron 
scattering and computer modelling study of MCM-41-S15 silica both dry and with 
absorbed light and heavy water, the latter samples giving a marked scattering 
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contrast due the change in neutron scattering amphtude between heavy and hght 
water. [301 Ell [22] According to Liu et al. [22] the diameter of the pores in this material 
is 15±lA. Two factors stand out from that study, namely, (a) the density profile of 
water in the silica pore is apparently not uniform across the pore, and (b) the density 
profile changes markedly when the temperature is lowered through 210K, being most 
uniform at that temperature and less uniform at temperatures either side of this value. 
The point here is that the measurement of a density minimum in confined water was 
achieved by measuring the intensity of the (100) Bragg peak in MCM-41-S15 silica as 
a function of temperature. |22] However the intensity of that peak will depend not only 
on the overall density of water in the pore but also on the way that water is distributed 
across the pore. Hence a change in intensity of the Bragg peak could equally be regarded 
as a due to a change in the overall density, or to a change in the distribution of density 
across the pore: this peak by itself cannot distinguish between these two scenarios. 

In this paper I show that the amount of water absorbed by MCM-41-S15, typically 
about 0.5gH2O/g substrate, is too large to be filling only the pores, assuming their 
nominal diameter of ISA. This means that either their actual diameter for absorbing 
water is larger than ISA or else some of the absorbed water appears on the surface of 
the silica particle and not in the pore itself. Either way, the discrepancy in pore sizes 
casts some doubt on the earlier observations of a fragile to strong dynamical transition 
in supercooled water, since water in the larger pore, or alternatively water on the surface 
of the silica particles, is likely to freeze on cooling, rather than go through a dynamical 
transition. Such surface water would of course, if present, be just as visible in scattering 
experiments as water in the pore. 

In addition I assess the intensity of the Bragg peaks in MCM-41-S15 silica as a 
function of water content (dry or fully hydrated) and of whether the water is protiated 
or deuteriated. As will be seen the availability of scattering data on each of these cases, 
namely dry, fully hydrated with H2O, and fully hydrated with D2O provides a significant 
constraint on the shape of the density profile. In particular it is shown that marked 
changes in the amplitude of these peaks, particularly for the deuterated material, will 
occur if the density profile of the water against the silica substrate changes, even if 
the overall mass of water present is unchanged. Hence the use of the amplitude of 
these peaks to assess the density of the confined water [22] may be unreliable. More 
importantly it is suggested that this analysis could be applied to other fluids absorbed on 
MCM-41-S15 silica and used to assess the likely surface density profile in those systems 
as well. 

2. Theory of Bragg scattering for a hexagonal lattice 

Figure [1] shows the first two atomic planes in the hexagonal lattice which give rise to 
the (100) and (110) Bragg reflections in Q (reciprocal) space. It is assumed here the 
cylinders extend indefinitely along their axes (perpendicular to the diagram) although in 
practice they are around l/xm long. Provided this linear extent is much larger than the 
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radius, as in the present case, the finite length of the real cylinder will not significantly 
affect the outcome of the present analysis. 




/ (100) / 

/ 

Figure 1. Schematic of hexagonal lattice of cylinders as found in MCM-41 silica. The 
lattice of hollow cylinders (shaded blue) is formed in an amorphous silica matrix. The 
axis of the cylinders is at right angles to the page. The first two Bragg planes, (100) 
and (110), are shown as pairs of dashed lines. The cylinders can be filled with various 
liquids, including water as in the present instance. This drawing is not necessarily to 
scale. 

Figure E] shows the actual diffraction data from the dry and wet MCM-41 silica and 
the Bragg peaks referred to above. If d is the spacing between cylinders, then spacing 
between the planes is dioo = 2 c? and duo = f • Since the reciprocal lattice vector for 
these reflections is perpendicular to the crystallographic c axis (which is parallel to the 
cylinder axis), these reflections can tell us nothing about the distribution of density along 
the cylinder axis. The Fourier transform of an inflnitely long, inflnitely thin cylinder 
will be a sheet of intensity in reciprocal space which lies normal to the axis cylinder. 
However the material in this case is polycrystalline with an assumed random orientation 
of crystallites, so the scattering intensity has to be averaged over these orientations. The 
Bragg peaks have to be further convoluted with the instrumental resolution function, 
which tends to be broad at low scattering angles. The Q values of these reflections will 
be -^ which works out, assuming d = 33. lA, to be 0.2191"^ for (100), 0.3791"^ for 
(110), and 0.438A~^ for (200). These values are close to the observed peak positions. 

In practice the cylinders are not inflnitely thin and will possibly have density 
inhomogeneities both along and perpendicular to the cylinder axis. Hence the Bragg 
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Figure 2. Differential scattering cross section for dry (a) and wet (b) MCM-41 silica 
showing the low Q Bragg peaks and the wider Q scattering used to estimate the atom- 
scale structure.^ Solid lines correspond to the 300K data, while the dashed lines 
correspond to 210K data. Note the logarithmic intensity scale. The inset in (b) shows 
the main water diffraction peak for the Wet-D sample at the same two temperatures. 
Note the marked shift of this peak to lower Q at 210K compared to 300K and the 
appearance of small Bragg-like features at 210K. 



intensities will be modified by the form factor arising from the finite size of, and density 
variation within, the cylinders. For simplicity we assume the scattering length density 
inhomogeneities are axially symmetric and so are a function of displacement from the 
cylinder axis (perpendicular radius, r^erp = r^ + r?/) and distance along the cylinder 
{z) only. Hence we can represent these inhomogeneities as p{x,y,z) = p{rpf,rp,z). The 
scattering amplitude of the cylinder in reciprocal space is then written: 



C{Q) = I dxdydzp{x, y, z) exp (iQ ■ r) 



(1) 



where r 



+ r. 



'-perp 



+ Tz. Q can also be split into perpendicular and 



parallel components: Q = Q 



perp 



Q^. Hence Q ■ r = Q 



perp' perp 



COS a + QzZ and 

dxdy = rperpdrperpda, where a is the angle between Qp^rp and Vperp- For the (100), 
(110) and (200) refiections, QzZ is zero by construction, so the integral over z can be 
performed directly and there can be no sensitivity in the Bragg peaks to the dependence 
of p {vperp, z) on z. Hcucc we replace p {rperp, z) by p (rperp), which represents the density 
as a function of perpendicular radius across the pore and averaged over the length of the 
pore. This is the radial density function calculated in the Empirical Potential Structure 
Refinement (EPSR) simulations [30], and will be referred to here as the density profile 
across the pore. 

Of course the Bragg intensity itself will be related to the value of \C{Q(^hki))\^-, and 
to calculate this intensity the density profile across the pore has be convoluted with the 
hexagonal distribution of pores and averaged over the orientations of the crystallites. 
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The simplest way to achieve this [20] is to set up a scattering length density function, 
nijc) = (6(r))p(r), where {b{r)) is the average scattering amplitude and p(r) is the local 
atomic number density at position r in the substrate or pore, and then perform the 
autocorrelation of n(r) to form a scattering cross section radial distribution function, 

G{r) = /[ n{r')n{r - r')dr\ , (2) 

where the angle brackets correspond to the average over the orientation of r. The 
scattering cross section is then calculated by Fourier transforming this to Q space: 

I(Q) = ^ f rG{r) sin Qrdr (3) 

This is the procedure adopted in the present work. For pure bulk silica (Si02), {b{r)) = 
5.25fm/atom, for H2O (6(r)) = — 0.56fm/atom, and for D2O {b{r)) = 6.38fm/atom. 
Assumed average atomic number densities will be 0.0663 atoms/A^ and 0.086 atoms/A^ 
for silica and water respectively. Silanol groups are attached to the inside of the 
pore, and these will be treated as OH groups {{b{r)) = 1.03fm/atom for OH and 
(6(r)) = 6.24fm/atom for OD). The fraction of Si atoms that have an OH attached 
is believed to be around 0.2 and is determined by inspection of the data as discussed in 
[30] and below. 



3. Experiment 

This paper will focus only on the intensities of the (100), (110) and (200) Bragg peaks 
for the wet and dry MCM-41-S15 using both H and D substituted water and for the 
two temperatures 300K and 210K. The experiment itself has been fully described in the 
preceding papers, [30l |3T1 132] and involves measuring the wide angle neutron diffraction 
pattern, over a wide Q range of 0.07 to 50A~^, from a series of dry and wet porous 
silicas, MCM-41-S15, whose nominal pore radius is 15A. The wet and dry samples 
are measured alternatively with H then deuterium (D) substitution on the water and 
silanol groups. In order to get the height of the Bragg peak onto a common scale for 
different water contents and different temperatures, the measured diffraction data are 
here renormalised onto a scale of differential cross section per scattering unit, with the 
scattering unit defined by the dry MCM-41 samples, namely csi = 0.28, cq = 0.61 and 
ch = 0.11, where Ca is the atomic fraction of component a. (See Table 1 in pQj.) The 
choice of scattering unit to normalise to is not particularly important for the present 
study. What is important is that it is the same scattering unit for all samples, H or D, 
wet or dry. That way the relative heights of the scattering levels and Bragg peaks can 
be compared between wet and dry samples. 

The renormalised data are shown in Fig. [2l Several features are notable from 
these graphs. Firstly, the high Q scattering level is determined by the composition of 
each sample. Since all the data have been normalised per atom of substrate, as the 
composition changes, so the high Q scattering level changes. There is a substantial 
change in the incoherent neutron scattering level between H and D isotopes (6.53 
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Table 1. Relative heights of the (100) Bragg peak for MCM-41-S15 for dry and wet 
hydrations and at 3G0K and 210K as a function of hydrogen isotope, H or D. All 
intensities are normalised to the intensity of the dry-D sample. 



Sample 


300K 


210K 


Dry-D 


1.00 


1.00 


Dry-H 


1.03 


1.03 


Wet-D 


0.21 


0.12 


Wet-H 


1.40 


1.32 



barns/sr/atom for H, 0.61 barns/sr/atom for D, 1 barn = 10~^^ni^), and as a result all 
the samples containing H have a significantly higher high Q scattering level compared 
to those with D. The H component in the dry samples arises from silanol groups on the 
pore surfaces. In fact the change in this scattering level between H and D isotopes in 
the dry samples is consistent with an assumed silanol fraction of ?^0.2, as discussed in 
[30] . The change in the scattering level between dry and wet samples is consistent with 
a water adsorption of approximately 0.3g water /g substrate for both H and D samples. 
This is somewhat lower than the 0.43g/g stated in [301 ES], but as will be seen below 
a discrepancy can arise because of the assumed radius of the pore, and the different 
numbers quoted here and in [30] may not be outside the experimental uncertainties in 
determining this value. 

A second observation is that in the dry samples, the main (100) Bragg peak barely 
changes in amplitude between the H and D samples, while the same peak in the wet 
samples changes drastically, with intensity in the D sample dropping by a factor of about 
7 compared to the corresponding H sample at the same hydration level. This large drop 
in the wet-D sample arises because the difference in scattering length density between 
pore region and silica region is large for the H sample, ~ — 0.3fm/A^, but smaller, 
•^ 0.2fm/A^, for the D sample (and the Bragg intensity will of course be proportional 
to the square of this amplitude). However the difference in scattering length density is 
not the only factor that contributes to the height of the Bragg peak, as will emmerge 
shortly. Table [T] lists the relative intensities of the (100) peak as a function of hydration, 
isotope and temperature, after substracting the underlying background from interfacial 
scattering, as discussed in [22] . 

At the same time the (110) and (200) peaks change with the different hydration 
levels and isotope substitutions. In particular the intensities of these peaks is markedly 
weaker for the D samples compared to H, particularly for the dry-D samples. (Note the 
logarithmic intensity scales in Figure [2] when comparing different data.) 

Hence these three aspects of the Bragg diffraction peaks as a function of hydration 
and isotope need to be reproduced by any model of the water density profile. 



Pw[r) 
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4. Data analysis 

4-1. Modelling the density profile 

For the purpose of the present discussion the water density profile is divided into three 
regions, namely a 'core' region, < r < re, representing water away from the silica 
interface, an 'interfacial' region, re < r < rj, and an 'overlap' region, rj < r < rp 
where water is interacting with the OH groups on the silica surface. The water density 
is specified separately for each of these regions, but is otherwise assumed uniform within 
in each region: 

Pc, for < r < re 

Pi, ioi re <r <rj 

(4) 
po, for rj<r<rp 

0, for rp < r 

The model therefore contains no atomistic information and cannot be used to calculate 
the structure functions at larger Q values. This choice of three water regions is intended 
as a simplistic representation of the actual density profile found in [30l [32] . The choice 
is certainly not unique, but is probably the minimum required to represent the actual 
density distribution obtained in [30] . 

The silica density profile is assumed to be zero in the core and interface regions, but 
can have a finite value in the overlap region, where it is treated as consisting entirely of 
OH groups (with the stoichiometry 0:H in the ratio of 1:2 to ensure charge neutrality). 
No Si is allowed in this region. This therefore assumes half of the oxygen atoms of 
OH groups do not occur in the overlap region but in the bulk region instead. The 
remainder of the volume is assumed to be the density and composition of pure bulk 
silica (8102). The density of the OH groups in the overlap region is obtained from the 
specified fraction of silanol groups (Si(0H)2) in the substrate as a whole, namely ~ 0.2, 
and with the assumption that no H atoms occur in the bulk silica region, and no Si 
atoms occur in the overlap region: 



Ps[r) 



(5) 



0, for < r < re 

0, for Tc < r < rj 

pOH, for rj <r <rp 

Pb, for rp < r 

where pp = 0.0663 atoms/A"^, as stated above. 

The above definitions also enable the respective volumes per unit length along the 
pore axis to be defined: 

(6) 

(7) 
(8) 
(9) 



Vc 


= TTr'c 




Vi 


= <r] - 


-rl) 


Vo 


= vr(4 


-r]) 


Vb 


= Vtotal 


— Trrp 
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where Vtotai = ^d'^ is the total volume of the hexagonal unit cell per unit length. 

If X is the silanol fraction then the average atomic fractions of Si, O and H in the 
substrate are given by csi = 1/(3 + 3x), cq = (2 + a;)/(3 + 3x), and ch = 2x/(3 + 3a:)[30]. 
If csi,B = 1/3 is the atomic fraction of Si in the bulk region and csi^o = is the assumed 
atomic fraction of Si in the overlap region, then the atomic density of the substrate in 
the overlap region can be calculated from: 

{cSi,B - Csi)pBVB ,^^. 

POH - —, TTT l-LUj 

[csi - csi,o)yo 

while the average density of the substrate is given by ps = iYoPoH + VbPb) / (Vo + Vb)- 
Note that setting the silanol fraction to zero in this model is equivalent to setting the 
volume of the overlap region to zero, since then csi^B = csi in f jTOj) and there is then 
no substrate density in the overlap region. In the analysis that follows the thickness of 
the overlap region, as determined by the difference rp — rj was treated as an adjustable 
parameter to give agreement Bragg peak intensities with experiment, while the silanol 
fraction is held constant at 0.2, in accord with the earlier experiment |30j . 

The outer radius of the overlap region, rp, is used to determine the total number 
of atoms of water in the pore, A^^^, = {pw)Vw, where V^, = Vc + Vi + Vo = Tfrj, per unit 
distance along the pore. 

In order to evaluate the radial density function, G{r) and hence the scattering 
intensity, I{Q), a hexagonal supercell is constructed consisting of 21 unit cells along 
each of the a and b crystallographic axes, with a cylindrical pore at the centre of each 
unit cell, and with an equivalent extension along the c axis which is at right angles 
to the hexagonal plane and parallel to the pore cylinder axes. Points are chosen at 
random within this supercell and assigned weights according to the substrate and/or 
water scattering densities at the position defined by the point, i.e. in the core, interfacial, 
overlap or bulk regions. The auto correlation function was then calculated directly from 
these weighted points out to a distance of r = ^lOrf, and, by using a standard Fourier 
transform, the corresponding scattering intensity in Q space is computed. Note that 
this simulation is not an atomistic empirical potential structure refinement (EPSR) 
simulation as performed by Mancinelli et al. [30] and so cannot reproduce details of the 
atomic and molecular structure. 

4-2. Estimating the pore radius 

The definitions of the previous section allow us to estimate the mass of water that can 
theoretically be adsorbed by different pore radii. If Aw = 6 is the average atomic mass 
per atom of water. As = J2a CaAa is the average atomic mass per atom of substrate, 
and we assume water occupies the pore uniformly, then the mass of water adsorbed per 
mass of substrate is given by 

M (rp) = {pw)7rrlAw/ (Vtotai - vrr^) pSAg. (11) 
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This function is plotted in Figure [3] for the case where the overlap thickness is held 
constant at 3.0A and the pore radius, rp, varied. Note that, as stated above, if the silanol 
fraction is set to zero (dashed line in Fig. |3]) this is equivalent to having zero substrate 
density in the overlap region. This figure demonstrates that neither the assumed overlap 
thickness nor silanol fraction makes a major difference to the conclusions about the 
possible pore radii for a given water mass adsorption. 



■I— < 
O 

(N 



1 



0.8 



0.6 



0.4 



0.2 







1 1 1 1 1 1 ' // 


t -a 1 .-. «-« *-» 1 4--«« .-. ^-» 4- -a *-» «-« / 1 ^1 1 


oiianoi iraciion u.z o / / 


t^ -< 1 r^ ♦^ J— * 1 \--ft-i jTi--« J— *«^ fill O , / 


oiianoi iraciion u.u / / 


Liu et al. ° o° // 

o / / 

// 
o / / 


o / / 
o /■■ / 

/ / 
/ / 


o / / 


o / / 


// 


O y^ / 


° -■ X 




o°° yy 


rP° 'V^ 


.y ^^^^ 


1 1 1 1 1 1 1 



7 



10 



11 
rp[A] 



12 



13 



14 



15 



Figure 3. Calculated mass of water per mass of substrate as a function of pore radius, 
for two silanol fractions, 0.0 (dashed) and 0.2 (solid). For the latter case it was assumed 
the inner radius of the overlap region, r/ extended 3 A below the pore radius, rp. Note 
that a mass adsorption of « 0.4g/g can only be achieved if the pore radius is of order 
12. 5A or greater. This calculation assumes the average water atomic number density 
is {pw) = 0.085atoms/A'^ as stated in [30]. Also shown is the case of Liu et al. [22] 
assuming zero silanol fraction and a pore separation of 25. 3A (circles). See text for 
more discussion. 



It will be noted immediately that if, according previous work [30j the mass of water 
adsorbed is of order 0.4g H20/g substrate, the radius of the pore must be of order 12. 5A 
(diameter ^ 25A) - given the pore spacing of 33.1 A there would not be enough volume 
for so much water in pores of significantly smaller radius. Moreover the composition of 
the substrate, i.e. whether silanols are assumed to be present or not, and whether or 
not there is an overlap region, does not affect this conclusion very significantly. 

One notable discrepancy between the data of Mancinelli et al. [30] and that of 
Liu et al. [22] is that for the latter case the (100) peak occurs at Q = 0.287A~^, not 
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Q = 0.219A^^ as here, even though the samples of MCM-41 are nominally the same. 
Such a large discrepancy in peak positions for the same material is surprising, but for 
the latter case would imply a pore separation of 25. 3A which is much smaller than the 
33.1 A pore separation in the present sample. Therefore the mass ratio for this smaller 
pore spacing is also shown in Fig. |3] as the dotted line. In that work the stated pore 
diameter is 15±lA meaning a pore radius of 7.5A, while the stated water adsorption 
is 0.5g D20/g substrate, which is equivalent to 0.45g H20/g substrate. However from 
Fig. [3] it is clear that such a large water absorption inside the pore is impossible for 
the specified pore radius of 7.5A, for which the maximum amount of water that can be 
absorbed is < 0.2gH2O/g substrate. A similar comment applies to other work, [5], where 
the pore diameter is stated to be 14A and the stated water absorption is 0.5gH2O/g 
substrate. 

There are two possible explanations of these discrepancies, namely either the pore 
radii are larger than is stated in any of the referenced papers, or else a significant amount 
of the absorbed water does not enter the cylindrical pores but is absorbed on the surface 
of the silica particles. Such surface absorbed water will most likely behave differently on 
cooling compared to the pore water, but will be equally visible in scattering experiments, 
and so can affect both the observed dynamic and Bragg diffraction behaviour. Hence 
it becomes rather important to establish the true size of the pores and the amount of 
surface water present. 

In fact there is another reason why, for the samples of Mancinelli et al. [30] , the pore 
radius may be of order 12. 5A. This is to do with the relative intensities of the (100), (110) 
and (200) Bragg peaks for the dry substrate. Figure H] shows the simulated scattering 
pattern for MCM-41-S15 with the pores empty for a series of pore radii encompassing 
the values given in Figure |3l Several factors compel us to insist that the correct pore 
radius for this material is near 12. 5 A. 

Firstly we can see that for smaller pore radii, the (100) peak in the dry-D sample is 
much weaker than for the dry-H sample: this is quite different from the behaviour seen 
in the measured data. Figure [2], where the (100) peak has almost the same amplitude for 
these two samples. Secondly in the data the intensity of (110) and (200) is significantly 
lower for the dry-D sample compared to the dry-H sample, with the (110) weaker than 
(200). This situation is only found in the simulation for a pore radius of 12. SA. Thirdly 
for pore radii smaller than 11. OA higher order reflections beyond (200) are found in the 
simulation for the dry-D sample which are not present in the diffraction data. Finally 
for a pore radius of 14. OA the (100) peak has markedly higher amplitudes for dry-D 
compared to dry-H. This trend is not found in the diffraction data, and the ratio of (100) 
to (110) and (200) intensities is much smaller than for smaller pore radii, suggesting 
again that a pore radius as large as 14.0 A is incorrect. 

The effect of changing the thickness of the overlap region, {rp — rj), using a fixed 
pore radius, rp = 12. sA and silanol fraction (0.2), is shown in Fig. |5l It is seen that in 
fact the sequence of Bragg peak intensities is rather insensitive to the overlap thickness 
at this pore diameter, although it will also be seen that for an overlap thickness less 
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Figure 4. Calculated Bragg intensity profiles as a function of pore radius {rp) , 
assuming a silanol fraction of 0.2 and an overlap layer thickness of 3.0A, for dry-H (a) 
and dry-D (b) MCM-41-S15 silica. The different pore radii are shown on the individual 
curves which are shifted vertically for clarity. 



than 3.0A, the (100) peak in the dry-D sample becomes progressively larger than the 
corresponding peak in the dry-H sample: this is the opposite of the experiment, where 
this peak is almost the same intensity in both samples. An equivalent graph (not 
shown here) was obtained when the pore radius was set to 8.0A, where the Bragg peak 
intensities were even more invariant to the overlap thickness. In that case however the 
ratio of (100) intensities for dry-H and dry-D samples is completely wrong compared to 
experiment as already shown in Fig. HI 

For a pore radius of 12. SA the effect of changing the silanol fraction at constant 
overlap thickness was negligible and so is not shown here, but for a pore radius of SA 
the effect is quite marked, as shown in Fig. O It will be noted that the (100) dry- 
H and dry-D Bragg peak intensities are only the same for this pore radius (as in the 
experiment) when the silanol fraction is set to zero. Since silanol groups are essential 
for water absorption in MCM-41 this seems like an unlikely scenario. 

The above analysis highlights the sensitivity of the Bragg peak intensities to the 
choice of pore radius, overlap thickness and silanol fraction. Combining the evidence 
of amount of water adsorbed per mass of substrate with this diffraction evidence, it 
seems pretty certain that the pore radius in the present sample of MCM-41-S15 was 
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Figure 5. Calculated Bragg intensity profiles as a function of overlap layer thickness 
{rp — rj) for a pore radius (rp) of 12. 5A, assuming a silanol fraction of 0.2, for dry-H 
(a) and dry-D (b) MCM-41-S15 silica. The overlap layer thicknesses are shown on the 
individual curves, which are shifted vertically for clarity. 



close to 12.5±0.5A. This is in contradiction to previously stated values for this radius. 
Traditionally estimating the pore radius in these porous sihcas has been a matter of 
some uncertainty. [33] . Almost certainly this arises from the variety of methods that 
are invoked to estimate this radius, which usually involve measuring the absorption 
isotherms of different gases adsorbed by the surface and using this information to 
estimate the surface area per unit length of the pore. This amount will depend 
intrinsically on the nature of the gas-surface or liquid-surface interaction, which can 
vary markedly from material to material. On the other hand the present diffraction 
evidence from the dry substrate, relies on no such assumptions, and so should be more 
reliable. The same diffraction method could be used to estimate the average pore radius 
in other porous materials where the pores are arranged in a regular lattice. 



4.3. Density profile in wet MCM4I 

Using only the diffraction evidence from the three Bragg peaks, (100), (110) and (200), 
probably it would therefore be difficult to guarantee a unique solution in determining the 
density profile. In this case however the earlier empirical potential structure refinement 
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Figure 6. Calculated Bragg intensity profiles as a function of silanol fraction x for 
a pore radius (rp) of 8 A and overlap thickness of 3 A for dry-H (a) and dry-D (b) 
MCM-41-S15 silica. The silanol fractions are shown on the individual curves, which 
are shifted vertically for clarity. 



(EPSR) simulations can be used as a guide [301 132] • In those simulations the water 
density at 300K in the interfacial region was observed to be ~2 times larger than the 
water density in the core region, with a decaying density in the overlap region. On 
the other hand at 21 OK the ratio of interfacial to core regions densities dropped to 
~1.5 under the constraint that the overall water density in the pore was not allowed to 
change. Are these variations consistent with the observed changes in the Bragg peaks 
between 300K and 210K? 

To answer this question the pore radius was set to rp = 12. 5A as in the simulations 
described above, with the core region extending to re = 6.5A and the interfacial region 
to r/ = 9.5A. The overlap thickness was therefore 3.0A. The silanol fraction was set 
to 0.2 as above, while the average water density was set to (pw) = 0.07 atoms/A^ for 
both 300K and 210K to give an overall water adsorption of 0.32g H20/g substrate. 
Although lower than the stated value in [30], this value is consistent with the observed 
change in high Q scattering level on going from dry to wet substrate (see Section [3]) 
and gives the correct simulated peak height dependence with temperature as observed 
in the diffraction data. At 300K the relative water densities were set to pc =1, pi =2, 
and po =0.7 for the core, interfacial and overlap regions respectively, with the actual 
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Table 2. Relative heights of the simulated (100) Bragg peak for MCM-41-S15 for dry 
and wet hydrations and at 3GGK and 21GK as a function of hydrogen isotope, H or 
D, calculated according to the simulations described in this paper. All intensities are 
normalised to the intensity of the dry-D sample at each temperature. 



Sample 


300K 


210K 


Dry-D 


1.00 


1.00 


Dry-H 


0.98 


0.99 


Wet-D 


0.20 


0.12 


Wet-H 


1.25 


1.24 



densities chosen to give the specified overall density. At 210K the same densities were 
set to 1, 1.5, and 0.9 respectively, giving a more uniform density profile than at 300K, 
as found in |32j. Figure [7] shows the results of these simulations and Table [2] gives the 
simulated relative peak heights to be compared with the measured values shown in Table 

m 

Clearly the qualitative trend of the diffraction data with temperature and between 
H and D samples is reproduced by these simulations. The main discrepancies are that 
the ratio of wet over dry intensity for the (100) peak in the simulation of the H sample 
is about 1.25, instead of 1.4 as observed in the experiment (see Tabled]), and for the wet 
D samples the heights of the (110) and (200) peaks are probably too large compared to 
the experiment. Hence some refinement of these values might be possible to improve 
the agreement with experiment. Figure [H] shows the subtrate density profile and the 
water density profiles at 300K and 210K used in this simulation. 

5. Discussion and conclusion 



At this stage an exhaustive study of all possible values of the density profile that might 
be consistent with the Bragg diffraction data has not been attempted. In addition 
only limited refinement of the values used in Section 14.31 has been performed, although 
such refinement should be straightforward. What does come out very clearly in this 
analysis however is that the height of the Bragg peaks, particularly those arising from the 
deuterated samples are highly sensitive to the choice of pore radius, overlap thickness, 
silanol fraction, and density profile across the pore, and so therefore cannot be used to 
state conclusively what is the density of water confined in porous silica, since a change 
in the density profile, such as that shown in Figure [8] can also have a marked effect on 
Bragg intensities. In particular the density profiles proposed in the earlier EPSR studies 
[30| 132] appear to be consistent with the observed Bragg peak intensities, within the 
approximations of the present work, which does not include an atomistic simulation as 
previously used. There is a difference between the present work and the previous work 
in that the assumed pore radius is here much larger (12. 5A) than assumed previously 
(9.6A). However in the earlier work there was significant penetration of the water into 
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Figure 7. Simulated Bragg intensity profiles for MCM-41-S15 silica assuming a silanol 
fraction of 0.2, for wet-H and wet-D materials at 300K (solid) and 210K (dashed). The 
H data are shifted vertically for clarity. Instrumental broadening, which is likely to 
be Q dependent and significant due to the low scattering angle at which the data 
have to be measured, as well grain size effects, single atom scattering, and pore size 
distribution effects, have not been included in these simulations. 



the silica matrix, which is covered by the overlap region in the present model, so in 
practice the two models may not be so different. Whatever the correct values for these 
parameters it is clear much more work needs to be done characterising the pore before 
any conclusions can be drawn about the nature of the water absorbed within it. 

Another feature of the literature on confined water is the (often) implicit assumption 
that the density profile in these systems is uniform across the pore. The highly charged 
nature of the silica surface, the fact that the surface is usually disordered, and the 
method of formation of these materials from cylindrical micelles, make this statement 
implausible at the atomistic level. Nonetheless evidence for this assertion comes from 
experiments in which the ratio of H20:D20 is chosen so that the scattering length density 
of the so-called contrast matched water is the same as the underlying silica matrix. 
Assuming the substrate is pure Si02 this ratio is approximately 0.4H2O:0.6D2O, and 
it is observed that within statistical accuracy the Bragg peaks disappear under these 
conditions [35]. It will be noted however in that work, that the (100) Bragg peak occurs 
at 0.17A~^, which would imply a pore separation of 42. 7A - significantly larger than any 
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Figure 8. Simulated atomic number density profile for dry MCM-41-S15 (bottom), 
for confined water at 300K (middle, shifted upwards by 0.1 atoms/ A"^) and confined 
water at 210K (top,, shifted upwards by 0.2 atoms/A'^). Vertical dashed lines demark 
the core, interfacial, overlap and bulk regions of the pore. 



values considered previously in this paper. In Fig. |9]the simulated Bragg intensities for 
the same contrast matched water using a pore separation of 42. 7A are compared with 
those from a simulation of wet-D MCM silica, using the same pore size (12. 5A), overlap 
thickness (SA) and silanol fraction (0.2) as in Fig. [71 It can be seen that although the 
simulated intensities do not go identically to zero, the Bragg peaks drop dramatically 
in amplitude and might be statistically too small to be observable. This occurs despite 
there being a highly non-uniform simulated density profile in the pore in this case, as 
in Fig. IHl Indeed the authors of this earlier work state that as water is desorbed from 
the pore, the density profile becomes non-uniform, with the highest density occurring 
near the edge of the pore, as in the present case. Since the previous experiments on 
water in confined geometry [30j were performed below saturation conditions, the finding 
of a non-uniform density profile in that case is perfectly reasonable. Here I have shown 
that negligible Bragg intensity with contrast matched water cannot be used to claim 
the density profile is uniform. 

A serious discrepancy with the earlier work of Liu et al. |22] is the pore radius, 
which as discussed in Section 14.21 is significantly different here to that stated in [22] for 
ostensibly the same material. The pore diameter determined from the present diffraction 
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Figure 9. Simulated Bragg intensity profiles for MCM-41 silica with the same pore 
separation (42. 7A) as given in [3S] assuming a pore radius of 12. SA, an overlap thickness 
of 3A, and silanol fraction of 0.2, for contrast matched water (bottom, H:D ratio 0.4:0.6) 
and D2O (top). The latter curve is shifted vertically for clarity. The density profiles 
across the pore and used in this calculation are taken from Fig. |8] at 300K. 



data, 25A, is near the nominal limit of pore diameters at which water will show a feature 
in the differential scanning calorimetry trace corresponding to freezing on supercooling, 
|34] . Nonetheless close inspection of the wet-D sample at 21 OK in the region of the 
main water diffraction peak near Q = 1.8A~^ (see Figure |2l inset) reveals three small 
peaks which appear closely analogous to those found in crystalline ice j3l], and other 
relatively sharp peaks appear at higher Q values in the 210K data which are not present 
in the 300K data. (See also Figure 2 of |32].) Indeed the shift of the main diffraction 
peak in D2O from Q ^ 2.0A-^ at 300K to Q ^ LSA^^ at 210K is similar to what is 
observed in heavy water when it freezes. Given the pore size broadening effect described 
by |3T], these peaks would be substantially broadened for a pore diameter of 25A. If 
the same crystal-like ordering transition occurred in the materials used to study the 
dynamic cross-over in supercooled water, then that cross-over would have a completely 
different explanation to the one normally advanced. This work of course does not rule 
out the possible dynamic crossover scenario to explain the apparent fragile to strong 
behaviour on cooling water absorbed in MCM-41 (or other substrate materials), but 
it does raise serious questions about the interpretation of existing experimental data 
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such as that hsted in the Introduction. The necessary pore size and wider Q atomistic 
information on the samples used to absorb the water is often unavailable: without that 
information the claims for observation of a FS transition or true density minimum are 
unsubstantiated, at least in the present case of MCM-41. Similar characterisation studies 
will needed for other substrates before similar conclusions can be drawn. Discrepancies 
to emerge in the present work on MCM-41 are that the pore separation clearly varies 
markedly from sample to sample, and there are significant doubts about what is the 
correct pore diameter: these parameters will have a profound effect on the amount of 
water that can be absorbed by the substrate. 

Quantitatively the density profiles presented here in Figure [H] are not identical to 
those obtained in the recent EPSR simulations, [32], but have considerable similarities. 
The main differences are in the overlap region where more water appears to penetrate the 
substrate in the present simulations than was found previously. In particular this overlap 
water increases on cooling to 210K whereas in the EPSR simulations it apparently 
decreases. Whether this difference can be resolved by formulating a version of EPSR 
that specifically fits the low Q Bragg peaks as well as the wider Q remains to be seen, 
but qualitatively at least it is correct to say the previously simulated density profiles 
are consistent with the observed change in Bragg peak intensity with hydration and 
hydrogen isotope. 

The present results of course do not categorically rule out the possibility of there 
being a density decrease on supercooling water in MCM41 silica. It is now clear however 
that the observed trends in Bragg peak intensities with temperature can be understood 
in terms either of a density change, or a change in density profile, or a combination 
of these two scenarios. What the present analysis demonstrates is that observation 
of a change in Bragg peak intensity with temperature does not, by itself, provide 
conclusive evidence for a change in density. It could equally be explained as a change 
to a more homogeneous density profile on lowering the temperature, and this latter 
scenario is consistent with previous atomistic simulations of the wider Q data at several 
temperature values. Hence the evidence for a density change is at best ambiguous at 
this point in time. A corollary of the density decreasing with cooling scenario is that 
if the density decreases, then the excess water will need to be excluded from the pore. 
This water can then freeze when it is eliminated, but will be perfectly visible in neutron 
scattering or other dynamics experiments. Such frozen water could also modify Bragg 
peak intensities, but more importantly might give a false signal of a dynamical transition 
when it is actually a freezing transition. 

More generally, for determining the scattering characteristics of both wet and dry 
ordered porous silicas, it appears that the present approach of simulating the Bragg 
peak intensity profile can provide invaluable information about both the pore size and 
the distribution of atomic density within it. This information is essential if we are to 
make accurate assessments of the nature of dynamic and structural transitions in the 
pore fluid. 
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